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ABSTRACT 


Equations are presented which govern the dynamics of the lines- 
first parachute unfurling process, including wave motion in the para- 
chute suspension lines. Techniques are developed for obtaining numeri- 
cal solutions to the governing equations. Histories of tension at 
the vehicle obtained vising the techniques are compared with flight 
test data, and generally good agreement is observed. Errors in 
computed results are attributed to several areas of uncertainty, the 
most significant being a poorly defined boundary condition on the wave 
motion at the vehicle-suspension line boundary. 
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INTRODUCTION 


The successful operation of a parachute depenas not only on a 
proper inflation phase, during which large aerodynamic drag forces 
are generated by the inflating parachute, but also on a proper phase 
of deployment, during which the parachute suspension system and 
canopy are extended downstream of the towing vehicle to a stage at 
which the inflation phase can begin in a satisfactory manner. In a 
general configuration, the parachute suspension lines and canopy are 
folded and packed into a deployment bag which is stowed in the towing 
vehicle until desired trajectory conditions are reached. The bag is 
then deployed rearward from the vehicle, either by some forced- 
}ec + ion device such as a mortar or by a small drogue or pilot para- 
lute. As the bag travels rearward, the packed suspension system and, 
thin, canopy sure unfurled from the bag. A typical vehicle-parachute 
configuration during a lines-first deployment is Bhown in figure 1 . 

Interest in such decelerator systems for planetary entry 
vehicles has stimulated research in deployment dynamics which has 
resulted in significant improvements over earlier empirical methods 
(references 1 and 2). A primary need for simulating the deployment 
process is to determine the motion of the deployment bag relative to 
the towing vehicle. Knowledge of sequence of events and unfurling 
times, the rates at which parachute material is unfurled from the 
deployment bag, and the levels of tension generated in the suspension 
system during the process of unfurling are important from mission 
analysis and deployment system design standpoints. 
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Reference 3 presents on extencive iiscussion of deployment 
systems and methods for analyzing deployment dynamics. Though res- 
tricted to the case of an inextensiblf. parachute, the techniques 
developed in reference 3 provide accurate simulation of vehicle 
trajectories and histories of deployment bag motion. The techniques 
also can be used to calculate reasonably accurate suspension-system 
tension histories, provided that averaged (as opposed to precise) 
parachute linear-mass-density profiles are employed. 

The results presented in reference U indicate that suspension- line 
elasticity is an important parameter in the prediction of tension loads 
experienced during parachute inflation. The significance of elasticity 
during the lines-first deployment process is indicated by the analysis 
presented in reference 5* Through a steady-state solution to the 
linearized one-dimensional wave equation governing suspension-line 
elasticity, a technique was developed in reference 5 which provides 
quite accurate prediction of the snatch force generated a c the paru- 
chute canopy begins to unfurl from the deployment bag. The conse- 
quences of neglecting wave mechanics in the modeling of suspension- 
line elasticity during the unfurling process were studied in reference 
6. It was found that, by using a massless-spring-type model for the 
elasticity, quite accurate histories of vehicle and deployment bag 
motion could be calculated, as was the case using the inextensible 
model. Generally accurate histories of suspension-system tension 
could be calculated except during periods of rapid load fluctuation 
caused by abrupt changes in the linear mess density of successive 




Figure i.- !>„ -.loyner.t configuration. 
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canopy mass elements unfurling from the deployment „;ig. Lai ge 
fluctuations in material unfurling rate corresponding to these linear- 
mass-density changes wer; calculated, but the accuracy of these cal- 
culations vas questioned due to the relative inaccuracy of the tension 
calculations. As a general conclusion, the inaccuracies were attributed 
to neglecting wave mechanics in the formulation of the model. 

The present paper employs the equations developed in reference t 
for governing the two-body planar motion of a lines-first forced- 
ejection deployment in conjunction with the one-dimensional nonlinear 
wave equation governing suspension-line elast response. In light o' 
the solution accuracy discussed in reference 6, the unfurling process 
is treated as two phases: a suspension-line unfurling phase, during 

which the formulation can be reduced to the massless-spring class; 
and a canopy unfurling phase, during which a solution to the complete 
set of equations must be obtained. Techniques are developed for 
obtaining a numerical solution to the governing equations. Generally 
good agreement between calculated results and flight test data is 
obtained by using the present model to simulate deployment loads and 
motions of two disk-gap-band par ac hut** deployment flight tests. 
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ANALYSIS 


In the analysir to * v ' 1 low, equations are presented which govern 
the dynamics of lines-first deployment of a parachute having elastic 
susp* nsion lines. Equations are presented which govern the motion of 
the towing vehicle and the motion of the deployment hag relative to the 
vehicle. An expression relating the tensio- developed at the mouth of 
the deployment bag to the rate at which the parachute exits the mouth 
of the bag (the unfurling rate) 1b given. The equation and boundary 
conditions governing the longitudinal wave motion in the elastic 
suspension lines are developed. In addition, expressions are developed 
by which the elastic behavior of the suspension lines can be approxi- 
mated by a massless-spring-class model during the phase of suspension- 
line unfurling. 

For the purposes of this analysis, the vehicle and the deployment 
bag and its instantaneous contrnts are consider?* to be mass particles. 
The surface of the planet is considered to be flat, and surface- 
relative accelerations are considered to be inertial The vehicle is 
assumed to be non-lifting, and its motion is restricted to a 
vertical plane. The motion of the deployment bag relative to the 
vehicle and the orientation of the tension vector in the parachute 
suspension lines are assumed to be parallel to the vehicle relative 
wind. 
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Equations of Motion 

The forces affecting the motion of the vehicle and deployment 
are snovn in Figure 2. 

The motion of the vehicle is influenced by aerodynamic drag, 
tension in the parachute suspension system, and gravitational attrac- 
tion. Under the previously mentioned assumptions, an expression 
was developed in reference 6 which governs the acceleration of the 
vehicle in the direction of its velocity vector, which is 


V 

V 


'VUt'v 


m 


+ g sin y 


( 1 ) 


In addition, the following trajectory equations are required in order 
to specify fully the earth-relative planar motion of the vehicle 


h “ v y sin y (2) 

;._g-£2 iX (3, 

v 

The motion of the deployment bag is influenced by aerodynamic 
drag, an unfurling resistance force, and gravitational attraction. 

An expression was developed, also in reference 6, which governs the 
acceleration of the deployment bag in the direction of the vehicle 
velocity vector, vhich is 


I I 
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n(C D A) b^ : F re 


+ g sin y 


(U) 


In this expression, the reduction in dynamic pressure acting on the 
bag due to the motion of the bag relative to the vehicle is also 
neglected. The velocity of the deployment bag relative to the 
vehicle can then be written, subsequent to solution of equations (l) 
and (h), as 


*b * % - v b <*> 

An additional expression relating the tension developed in the 
unfurling parachute at the mouth of the deployment bag to the linear 
mass density of the parachute element exiting the bag and the unfurl- 
ing rate was given in reference 5» and can be written as 

T_ = m'u 2 + F (6) 

B re 

The unfurling rate is, by definition, the rate of change of the 
spatial coordinate of the mouth of the deployment bag, or 


u ** 


(T) 
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Suspension-System Elasticity 

Complete definition of the motions of the deployment process 
requires the formulation cf some relationship between the tension 
generated at the mouth of the deployment bag , the tension transmitted 
to the vehicle, the relative motion of the deployment bag, and the 
unfurling rate. Such a relationship is furnished by a mathematical 
model of the parachute elasticity. In the present paper, the parachute 
canopy is considered to be an inextensible structure due to its complex 
cloth-tape configuration, which is thought to undergo very little, if 
any, deformation during t.he unfurling process. The dynamic response 
of the suspension lines during the unfurling process can be determined, 
most accurately by a mathematical model which considers the wave 
mechanics of the lines (reference 6). As the deployment motion is 
considered 'n the present paper to be one-dimensional, such a model 
requires only a single partial differential equation to govern the 

•X/ 

longitudinal wave motion in the lines. 

The present model considers the parachute suspension system to 
consist of suspension lines only, which have negligible aerodynamic 
drug and uniform elastic and mass properties. The elastic 3tate of 

the suspension lines can be represented by use of the coordinates t 

and variables shown in Figure 3* The governing equation for the * 

longitudinal motion of an arbitrary suspension-line element was • 

developed in reference 6 and can be written, in light of the present 

\ 

"as mptions, as 
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3T 
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+ T„ 
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The elastic behavior of the element can be described by a relationship 
between the tension, T, and the displacement, 6, and its derivatives. 
Considering the tension to be a sum of static and dynamic components, 
this relationship can be written as 


■ V (> & * e St (If)] 


( 9 ) 


The present model considers the specific secant modulus, K' , to "be 

sec 

36 

a function of the strain, |y , and assumes the damping coefficient, C, 
to be constant. 

Boundary conditions at the point at which the suspension lines 
attach to the vehicle and at the mouth of the deployment bag (or at 
the Juncture of the suspension lines and the canopy, for times after 
which the lines have been completely unfurled) must also be prescribed. 
The present paper considers the suspension lines to be attached to the 
vehicle in a manner such that no deflection in the lines is allowed 
at that point. Mathematically, this condition can he stated as 






u. 


do) 
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The tension is subject to the boundary condition at the mouth of the 
deployment bag which is given in equation (6). Additionally, an 
equation of compatability between the motion of the deployment bag 
and the motion of the exiting parachute element can be written as 


D6i 

Dt 1 


1=1 


B 


u, 


(0 




L sl> 


( 11 ) 


Suspension-Line Unfurling Phare 
It was concluded in reference 6 that a massless-spxxng-type 
model of the suspension-line elasticity provides acceptable simulation 
of deployment loads and motions during periods of small load 
fluctuations, as is the case during the suspension-line unfurling 
phase. In order to avoid the numerical complexities of solving a 
partial differential equation in an expanding space, a similar model 
is employed in the present analysis to govern the suspension-line 
unfurling phase. Under the assumption of a massless-spring model, 
equation (8) reduces to 


ft ' 0 • < 12 > 

Therefore the tension is constant throughout the lines at any 
particular time, or 


T„ * T. 


B 




sr 


(13) 
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The total derivative in equation (11) can he expanded into local 
and convective terms and, considering the definition of the unfurling 
rate as given hy equation (7), can be written as 


D6 i 96 i 96 i 

Dt = 0+- 1 u 1 

«,=«. dt dXr 1=1 

B B B 


(lM 


Since the strain rates experienced during suspension-line unfurling 

are very small, the dynamic component of the tension expressed by 

equation (9) can be neglected. Then equation (9) can be solved for 

36 

the strain , which can be evaluated at &=)lg and substituted 
into equation (l4)» resulting in 


D6 = 96 | 

Dt 9t ^ . 

13 


+ u 


n ,K' 
si sec 


(15) 


Tp can be replaced in equation (15) by its definition as given in 
equation (6) and the entire expression substituted into equation (11). 
This results in a cubic expression for the equilibrium unfurling 
rate at any time, which is 



36 

However the deflection rate at the bag mouth, rr-| , is several orders 

U£ B 

of magnitude smaller than the relative bag velocity during the 
suspension-line unfurling phase. Thus, equation (16) can be reduced 


y/ 


-> 



to 



1U 



Equations (l) through (T). in conjunction with equations (13) and (17), 
thus govern the suspension-line unfurling phase. 

Canopy Unfurling Phase 

The canopy unfurling phase begins when the suspension lines are 
completely unfurled from the bag. As rapid load fluctuations 
characteristically occir during this phase, a solution to the equation 
governing wave motion in the suspension lines must be obtained (equa- 
tion (8)). Initial conditions on the motion of the suspension lines 
at the beginning of this phase are prescribed by the massless-spring- 
type solution obtained during the previous phase. As the tension 
(and thus strain) at any particular time is uniform throughout the 
length of the lines in that solution, the initial displacement profile 
must be linear (with = 0) in order to assure coinpat ability 

between the phases. The initial displacement rate profile is assumed 
to be identically zero, in accordance with the assumption in effect 
during the previous phase. During the canopy unfurling phase, the 
boundary condition expressed by equation (11) is applied at the 
spatial coordinate of the Juncture of the suspension lines and canopy, 
which is stationary. Thus, the boundary condition can be reduced 
to partial derivative form and can be expressed as 



15 
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3t 


£.=L 


h - u > (L el < l B i V 


(16) 


A simultaneous solution of equations (l) through (9). in conjunction 
with the boundary conditions given in equations (10) and (l8) and the 
initial conditions described previously, defines the loads and 
motions of the canopy unfurling phase. 



TECHNIQUES FOR NUMERICAL SOLUTION OF GOVERNING EQUATIONS 

In order to define the loads and Lotions of both the suspension- 
line and canopy unfurling phases, a solution to the ordinary 
differential equations governing the motion of the vehicle and deploy- 
ment bag must be obtained. During the canopy unfurling phase, a 
solution to the partial differential equation governing wave motion 
in the suspension lines, subject to prescribed initial and boundary 
conditions, must also be obtained. Techniques for obtaining proper 
solutions during both unfurling phases have been developed using 
standard finite-difference methods, such as those described in 
reference 7, and are described in this section. 

Suspension-Line Unfurling Phase 

The suspension-line unfurling phase is governed by six first-order 
ordinary differential equations (equations (1) through (5) nnrl ( 7 ^ 
in conjunction with an algebraic equation for calculating the unfur] - 
ing rate (equation (17)) and two expressions for calculating the 
tension generated at the mouth of the deployment bag and the tension 
transmitted to the vehicle (equations (6) and (13)). Initial values 
of the variables ( , h, y, v^, x^, and Jig) of couree depend on 
flight trajectory conditions and the individual system configuration. 

The set of differential equations is integrated using the single- 3tep 
forward-difference technique suggested by Euler. Given the value of 
the first derivative of a time- variant function y at time station i, 
Euler's technique is used to calculate y at time station i+1 as follows 


1 6 
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y i+i * y i + (19) 

This procedure is quite simple and is particularly useful for integrat- 
ing smooth functions such as those derivatives of interest during 
this phase. Since the frictional resistance force F and the ri.cpen- 
sion-line linear mass density m^ are constants and the specific 
secant modulus is a smoothly varying function of the strain 

(and, actually, is nearly constant during this phase) equations ( 17 ) 
and (6) can be solved without iteration for the unfurling rate and 
tension at the deployment bag mouth, respectively. 

Canopy Unfurling Phase 

The motions of the vehicle and deployment bag during the canopy 
unfurling phase are again governed by the six first-order ordinary 
differential equations as described in the previous section. Values of 
the variables of concern at the initiation of this phase ore, of 
course, their values at the end of the suspeneion-line unfurling 
phase. These differential equations are again integrated using Euler's 
technique. 

A solution must be obtained simultaneously to the second-order 
partial differential equation governing wave motion in the suspension 
lines (equation (8)) in conjunction with the elastic model expressed by 
equation (9)* The wave equation is subject to the boundary condition 
of no deflection at the vehicle attachment point, as expressed by 
equation (10), and to the boundary condition of compatability at the 



Juncture of the suspension lines and canopy, as expressed by 

equation (18). Initial conditions on the wave motion are those 

discussed in the analysis section, namely, the initial strain profile 

is linear, with Mjj-q s 0» and 'the initial deflection rate profile is 

identically zero. The wave equation car. be solved by using stanaaru 

central-difference formulations which can be derived using the gric 

pe.tterr and subscripts shown in Figure U. 

The derivative — ;• at any interior point can be appro, imat'-d t.y 
3t‘* 

the following central-difference equation: 

A\ . Vuj (toi 

H i.j (it> 

3T 

The tension gradient at any interior point can be approximated 

by a central difference equation by using the half-station concept, 
as follows: 


9T 

di 


i.J 


T l-Jtl/2 

Li 


( 21 ) 


The tensions at the half stations i,J+l/2 and i,J-l/2 are a 
function of the local strain and strain rate, as expressed by 
equation (9)- The strain at stations i,J+l/2 and i,J-l/2 can 
be evaluated by using the following difference approximations: 


111 

H i,j+l/2 


^ j a , ^ i 

, Li 1 

LI 


( 22 ) 
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36, 

35L 1 


i»J-3/2 




( 23 .: 


The strain rates at these two stations can be evaluated by using the 
following difference approximations in conjunction with equations 
(22) and (23). 


3 _ ,36 
3t K dl 


i»J + l/2 


— (— )| 
3t K n n 


i,J-l/2 


36. _ 36. 

i ,J+l/2 dt i-l,J+l/ ~ 

At 


36.i _ 36. 

A i.J-1/2 31 i-l_, j-1/2 

At 


(2t 


( 25 ) 


A value of the specific secant modulus at each of the half-stations 
can be found explicitly as a function of the local strain. Then 
values of local strain, strain rate, and specific secant modulus 
car. be used in equation (9) to calculate the tension at the half- 
ntationa. Thus, the approximations given by equations (20) and (21) 
can be substituted into equation (8), resulting in a difference 
approximation to the wave equation which can be used for calculating 
future values of the deflection 6 at station j . This 
approximation is 


* 26 i.j - * l4t) ‘ 


n . 1 '. 
si si 


T 1 ..1+1/2 ~ T i.. 1-1/2 
A* 



(26) 
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In using the difference approximation given by equation (26), care 
must be taken in the selection of values of At and A£ in order to 
assure a stable solution. Stability can be assured if the ratio of 
AH to At is greater than the maximum wave propagation velocity 
anticipated in the course of the solution. Using the classical 
definition of wave velocity, the stability criterion which must be 
satisfied can be stated as 



where K' is the maximum value of the specific secant modulus 

sec, max 

in the range of suspension-line strains expected in the course of the 
solution. 

The deflection at spatial station n, which is located at the 
juncture of the suspension lines end canopy (8. = L g ^),m"st be obtained 
through a simultaneous solution to the wave equation (equation ) ) 
and the boundary condition of compatability expressed by equation 
(l8). In the development of a technique for obtaining this solution, 
it is assumed that the tension level is constant at all points through- 
out the canopy (acceleration of the unfurled portion of the canopy is 
neglected). Then the tension at the Juncture of the suspension lines 
and canopy equals the tension being geneiated at the mouth of the 
deployment bag, which is expressed by equation (6). It is also assumed 
that the linear-mass-density profile of the canopy can be represented 


§ 




mathematically in a manner such that no discontinuities are present. 
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The inter-relationship of equations (6), (8), ai.d (l8) suggests 
some type of iterative solution technique. Such a technique can be 
developed by first treating the wave equation, equation (8), as two 
first-order equations. By defining 


V 

n 


36 

3t 




si 


equation (8) can be written as 


(28) 



n .m' 
si si 
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si 


( C_A ) a + T 
D V 
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(29) 


The tension gradient in equation (29) can be approximated by the 
following backward-difference equation (again using the subscripts 
shown in Figure U': 


3T 

3£ 


i»n 



IkjL. -ALL 



( iO ) 


In this equation, T 0 is prescribed by equation (6), and T i,n-l/2 can 
be calculated using the difference procedures outlined previously. A 
first estimate of the integral of equation (29) can then be calculated 
using Euler's method as gx/en by equation (19) » or 


V 

n 


i+1 


V 

n 


+ At 
i 



(31) 
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Provided that an estimate of at time station i+1 is known, a 

first estimate of the unfurling rate at that time station can be 
calculated using equation (18), or 


u 


i+1 



i+1 


( 32 ) 


This value for the unfur] ing rate is used in equation (6) to calculate 
a first estimate of the tension at the hag mouth (and, t'us, at 
spatial station n) at time station i+1. However, since the linear 
mass density of successively unfurling elements can be varying quite 
rapidly, this value of the tension at the bag can be substantially 
different from the value at the previous time station. Therefore, it 
is necessary to iterate on the solution. Ideally, the iteration 
procedure should include solution of equation (26) (eo all spatial 
stations other than station n) and integration of equation (28), in 
order to update the array of deflection values 6 for a updated 
estimate of the tension at station n - 1. In consideration of the 
fact that such a series of calculations would require extensive 
computer time, and also that, for sufficiently small time steps, 
changes in the deflection array will be small in comparison to 
changes in , the deflection array is not updated until the iteration 
on the boundary solution has converged. 

In light of these considerations, the iteration procedure is 




structured as follows: 
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(l) The first estimate of T | , which has been calculated 

b i+1 

previously, is substituted back into equation (30) and an estimate of 
the tension gradient at time station i+1 is calculated as follows: 


9T. 

34 ' 


i+l,n 


T - T 

B'i+1 i,n-l/2 

1 “ 


(33) 


(2) This value of the tension gradient can be substituted into 

3v 

equation (29) and a value for can be calculated. If this 

3V q dX i+1 

new value tt— | is within some error bound of the previous value 
3V dt i+1 

^— | , the iteration procedure has converged, 
i 

(3) If the error criterio 1 ' is „ot satisfied, a new estimate of 
the integral of equation ( 28 ) can be obtained using the trapezoidal 
method of numerical integration, which can be written as 


i+1 




3V 


ni 


3t 


av 

E 

at 


i+1 


(3U) 


This new estimate of V | can then be substituted into equation 


i+1 


(32) to calculate a second estimate of the unfurling rate u ^ +1 * A 


second estimate of T. 


B' 


can then be calculated, again using 


i+1 


equation (6). The procedure is then repeated until two successive 

av 


values of 


at 


lie within the error bound of each other. 


i+1 


Since a converged value of V 


i+1 


has also been found, it can 


be used in a trapezoidal integration scheme to find as follows: 



W * * f 


+ V 


n' 


i+1 


(35) 
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The final parameter to be evaluated, is the tension at the point 
at which the suspension lines attach to the vehicle, T^. The deflec- 
tion 6 at this point is assumed to be equal to zero, as expressed by 
equation (10). The basic partial differential equation (equation (8)) 
also holds at this point. By ^ mbining these two equations, an 
expression for calculating the tension gradient at the vehicle attach- 
ment point can be derived. 

Since the deflection at the attachment point is always zero, its 

derivatives with respect to time, namely — t \ , are also zero. 

3t 4-0 

Thus, for this point, equation (8) can be reduced to 


22 a 

34 1 


= n ,m 


4=0 


si si 


<°d a) A- 


+ T, 


m 


( 36 ) 


An additional expression for this tension gradient can be found by 
exiianding the tension at the mid-point of the grid space between 
spatial stations 1 and 2 (station 3/2) into a Taylor's Series about 
station 1 (at which 4=0). Neglecting second- and higher-order 
terms, this expansion is 


T i,3/2 


A i,l 2 34' 


i.l 


(37) 


The tension at station i, 1,2 can be calculated numerically by using 
the difference procedures outlined previously, with Ihe subscript 
J = 1. Since the tension gradient in equation (37) c *n be evaluated 
directly by using equation (36), the tension at the attachment point, 
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T 1 , can then be evaluated directly as 


m - rn 

i,l 1,3/2 


AA 

2 n 8l iu sl 


tyyi + T v 


m 


(38) 


Ideally, since (which, numerically, is also appears 

on the right-hand side of equation (38), a scheme for iterating on 

this equation in some manner until and T\ ^ converge to the 

same value would be preferred. However, since m' , is generally an 

si 

extremely small number, appreciable changes in result in only 

minute changes in the magnitude of the tension gradient (equation (36)), 
and thus would alte. only slightly the value of ^ in successive 
iterations. Thus, equation (38) is solved only once at each time 
step in the present paper. 

The sequence of numerical calculations described in this section 


is outlined in Figure 5- 
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Figure 5.- Sequence of numerical procedures 
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RESULTS AND DISCUSSION 

A precise treatment of suspension-line wave motion during the 
parachute unfurling process is virtually impossible due to non- 
uniform mass and elastic characteristics, multi-dimensional motion, 
and poorly defined boundary conditions to which the motion in subject. 
The present method combines the planar representation of vehicle and 
bag motions with an idealized treatment of the wave motion in an effort 
to achieve reasonably accurate simulation of motions and loads In 
order to evaluate the accuracy of the present method, it was used to 
calculate deployment loads and motions for two flight tests of disk- 
gap-band parachutes: the second balloon-launched flight test of the 

NASA Planetary Entry Parachute Program, B/L-2 (reference 8) and the 
flight test of vehicle AV-U of the NASA Balloon Launched De^elerator 
Test Program (reference 9)- Physical system data and parachute mass 
distributions for the two flight tests are presented and discussed 
in the appendix. As there is an absence of data, a numerical value 
for the suspension-line damping coefficient was determined parametri- 
cally. As the suspension system never became slack during the two 
flight test deployments the damping coefficient selected was the 
smallest value required to numerically maintain tension at the vehicle 
attachment point during the unfurling process. 

Computed histories of unfurled length for the two flight tests 
are shown in Figure 6 and compared with test data points corresponding 
to mortar fire, line stretch, and estimated bag strip. The computed 
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histories exhibit very good agreement with flight test data. 

A sequence of computed profiles of tension in the parachute 
suspension lines is presented for the B/L-2 flight test in Figures 7, 

8, and 9- Figure 7a shows the state of the suspension lines just 
prior to the emergence of the skirt of th«* Tjarac’ ute canopy from the 
bag. The tension is constant along the length of the lines as is 
consistent with the assumptions made in the so lotion during the 
suspension-line unfurling phase. Unfurling of the lower edge of tne 
band from the bag generates a sharp peak in the tension profile, as 
is shown in Figure Tb. At, this peak travels down the suspension 
lines, it is damped and spread; Figure 7c shows the wave at> it strikes 
the vehicle. The reflected wave is shown traveling back up che 
suspension lines in Figure 8a. The unfurling of the top of the band 
section of the canopy from the bag generates a second peak in the 
tension profile, which is shown in Figure 8b. This wave again is 
damped ar.d spread as it travels toward the vehicle; the tension profile 
as the second wave strikes the vehicle is shown in Figure 8c . Un- 
furling of the bottom edge of the disk section of the • anopy from the 
bag generates a third wave, which is shown departing the Juncture of 
the suspension lines and canopy in Figure 9a* Figure 9t shows the third 
wave striking the vehicle as the reflected second wave reaches the 
Juncture. Although interaction of generated and reflected waves 
causes some fluctuation in the tension at the vehicle attachment 
point after the third wave has reflected, the tensior at that point 
decays rather smoothly due to the decreased tension being generated at 



Unfurled length, meters Unfurled length, meters 
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Figure 6.- Comparison of computed histories of unfurled length with 
flight test data. 
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Spatial coordinate, f f meters 
(c) First tension wave strikes vehicle (t^ • .659 sec.) 


Figure 7'~ Computed profiles of tension in suspension lines during 
FEPP B/L-2 parachute deployment. 
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Spatial coordinate, { , meiers 
(c) Second tension wave strikes vehicle (t^ ■ .732) 

Figure 8.- Computed profiles of tension in suspension lines during 
PEPP B/L-2 parachute deployment (t . * . 678 , .709, and 
.732 sec). 
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the bag mouth as the disk unfurls. A fourth wave (of smaller magnitude 
is generated upon the unfurling of the vent edge from the bag as is 
shown in the profile in Figure 9c. However, bag strip occurs before 
this fourth wave reaches the vehicle. 

Tension profiles very similar in nature to those shown in 
Figures 7, 8, and 9 were obtained in the simulation of the AV-U flight 
test. In the interest of brevity, only three of these are shown. 

Figure 10a shows the wave which is generated by the unfurling of the 
lower edge of the band as it leaves the juncture of the suspension 
lines and canopy. Figure 10b shows the second wave, which is generated 
by the unfurling of the top edge of the band, as it leaves the juncture 
and also shows the damped, reflected first wave traveling back toward 
the juncture. The profile shown in Figure 10c shows the third wave, 
which is generated by the unfurling of the lower edge of the disk, 
leaving the juncture as the second wave iB approaching the vehicle. 

Histories of the computed tension at the vehicle attachment point 
are compared with flight test data histories in Figure 11, for B/L-2, 
and in Figure 12, for AV-4. Peak loads shown in the initial half of 
each flight test history are attributed to phenomena which are not 
associated with the dynamics of unfurling and, as such, are not 
included in the present mathematical model. There is very good 
agreement in occurrence times of peak loads between the flight-test- 
data histories and the histories ccmputed using the present method. 
Although agreement in peak load magnitudes is only fair, the results 
computed using the present method show considerable improvement over 
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the results obtained using the massless-spring model (reference 6). 

Errors in the present computed histories can be attributed to several 
areas of uncertainty: 

(1) Lack, of data concerning the elastic response of suspension- 
line material at very high loading frequencies. 

(2) Probable damping of waves generated by the unfurling of upper 
sections of the canopy as they are propagated through 
sections of the canopy which have been unfurled previously. 

(3) Dispersions of wave propagation velocities between individual 
suspension lines. 

(U) A pocriy- defined boundary condition at the vehicle boundary. 

In flight tests, the suspension-lines are, in general, 
attached to a multi-legged, webbed bridle system, which in 
turn is attached in some manner to the vehicle itself. The 
bridle webbing probably exhibits rather severe damping 

\/ 

effects and, in addition, waves reaching the vehicle itself 
are probably partially absorbed rather than totally reflected. 

Computed histories of relative deployment bag velocity and 
unfurling rate for the two flight tests are presented in Figure 13. 

Similar trends are exhibited in both flights. The relative bag \ 

velocity decays smoothly due to deceleration of the towing vehicle. 

The general trend of the unfurling rate is a similar decay. • 

However, sharp decreases in the unfurling rate occur when sudden f 

increases in the linear mass density of the unfurling parachute 

s 

are encountered; conversely, when sudden decreases in the linear 
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Figure 12.- Comparison of computed and flight-test histories of 

tension at vehicle during BLDT AV-U parachute deployment. 
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mass density are encountered, the unfurling rate increases sharply. 

In comparing the computed B/L-2 unfurling rate history with the results 
presented in reference 6, it is evident that the magnitudes of the 
maxima and minima axe decreased and increased, respectively, in the 
present history. Specifically, the unfurling rate as the gap portion 
of the canopy emerges from the deployment hag exceeds the local 
relative hag velocity hut does not exceed the B/L-2 flight-test mortar 
ejection velocity, as the results in reference 6 showed. Although 
there is an absence of data concerning unfurling rate, the magnitudes 
presented in the present paper axe thought to he good approximations 
to the magnitudes experienced during the flight-test deployments 
being studied. 
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CONCLUSIONS 

The present paper has presented the equations governing the lir.es- 
first parachute unfurling process, including the equation and boundary 
conditions governing idealized longitudinal wave motion in the para- 
chute suspension lines. Numerical techniques for obtaining a solution 
to the governing equations have been developed. Based on numerical 
solution of the governing equations for two disk-gap-band parachute 
deployment flight tests, the following conclusions can be made in 
evaluation of the present methods of simulation: 

1. Computed histories of unfurled length agree very well with 
flight test data. 

2. Computed histories of tension at the vehicle agreed well, in 
general, with flight test data histories. Errors in the 
computed histories were attributed to several areas of uncer- 
tainty, the most significant of which is probably a poorly 
defined boundary condition on the wave motion at the vehicle 
boundary. 

3. Computed histories of unfurling rate exhibited fluctuation 
about the smooth histories of relative deployment bag 
velocity. However, the amplitudes of this fluctuation were 
attenuated as compared to the fluctuation experienced with 
the use of a massless-spring-type model of suspension-line 
elasticity. Specifically, no local values of unfurling rate 
computed using the present methods were observed to exceed 
the respective flight-test mortar ejection velocities. The 
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present computed histories are thought to be good approxima- 
tions to the unfurling rates experienced during the flight 
tests being studied. 
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APPENDIX 


PEPP B/L-2 AND BLDT AV-1* SYSTEM DESCRIPTIONS 
AND FLIGHT TEST DEPLOYMENT CONDITIONS 

Proper simulation of the PEPP B/L-2 and BLDT AV-U flight t«*st 
deployment requires, of course, correct numerical input data. These 
data include mortar-fire trajectory conditions ; vehicle mass, diameter, 
and aerodynamic characteristics; deployment hag mass, diameter, urag 
coefficient, and the velocity with which it is ejected by the mortar; 
parachute sizes and mass distributions; and suspension-line elastic 
characteristics. The data which were used to obtain the results 
presented in the present paper are described in thi6 appendix. 

Mortar-Fire Trajectory Conditions 
Trajectory conditions at mortar fire for simulation of the two 
flight tests were selected so as to match the most important (with 
respect to the dynamics of the deployment process) variables, namely, 
the vehicle velocity v y , the flight-path angle y, and the freestream 
dynamic pressure q^. Values of these initial parameters are listed 
in references 8 and 9- The present simulation technique utilizes 
tables of atmospheric data from the U. S. Standard Atmosphere, 1962 
(reference 10 ). Atmospheric data obtained prior to each of the flight 
tests were noticed to deviate slightly from those of the Standard 
Atmosphere. These deviations were taken into account by adjusting 
the input values for altitude at mortar fire downward to values at 
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vhJch the tabular densities matched those local data values obtains 
prior to the flight tests. Mortar-fire trajectory parameters which 
were used in the simulation are listed in Table I. 



TABLE I 


INPUT TRAJECTORY CONDITIONS AT MORTAR-FIRE 



PEPP B/L-2 

BLDT AV-l 

Altitude, h, kilometers 

39-8 

1U.1 

Vehicle velocity, v^, me". -rs/ second 

518 

698 

2 

Freestream dynamic pressure, nevtons/meter 

555 

522 

Mach number 

1.63 

2.13 

Flight-path angle, degrees 

50 

12 


Vehicle Descriptions 

The B/L-2 flight test vehicle consisted of a 120°-inciuded-angle 
cone aeroshell which enclosed a cylindrical payload. Mass, diameter, 
and drag characteristics (as a function of Mach number) for thi3 
vehicle were obtained from reference 11. The AV-4 flight test vehicle 
was a full-scale model of the Viking '75 entry vehicle. Mas3, 
diameter, and drag characteristics for this vehicle were obtained from 
reference 9* As the flight test Mach numbers for B/L-2 and AV-l 
varied only slightly during the unfurling process, the respective 
drag coefficients can be approximated by constants (see Table II). 
Tables of the wake parameter n as a function of distance aft of the 
vehicle bases were calculated using flow properties (namely, dynamic 
press’ire and flow velocity) at the center of the vake of the two 
vehicles, which are presented in reference 12. The wake parameter 
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varies from a value of 0 at the base of the vehicle to values near 
0.75 in the region of the wakes in which bag strip occurs. Vehicle 
masses, diameters, and approximate drag coefficients are summarized in 
Table II. 


f 




TABLE II 


VEHICLE PHYSICAL CHARACTERISTICS AND DRAG COEFFICIENT’S 

PEPP B/L-2 BLDT AV-k 

Reference diameter, meters L.6 3-5 

Mass, m^, kilograms 566 8 IT 

Approximate drag coefficient, C n 1.5 1.6 

D ,v 

Deployment-Bag Chfiracteristics and 
Mortar Ejection Velocities 

The parachutes for the respective flight t°sts were pressure- 
packed into cylindrical fabric deployment bags. Approximate bag masses 
were obtained from references 8 and 9» The bag diameter for each 
flight test was taken to be 0.3 meters. Drag coefficients for the bags 
were assumed to be constant and equal to 0.8, which is the approximate 
drag coefficient of a cylinder in uniform axial flow (reference 13 ). 

The unfurling resistance force F^ e was determined in laboratory 
tests to be equal to an essentially constant value during the suspension- 
line unfurling phase and another essentially constant value during the 
canopy unfurling phase. Mortar ejection velocities were determined 
parametrically as those values required to obtain agreement between 
computed and flight-test data times from mortar fire to line stretch. 
These values, along with the previously mentioned deployment bag 
characteristics, are summarized in Table III. 
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TABLE III 

DEPLOYMENT BAG CHARACTERISTICS 
AND MORTAR EJECTION VELOCITIES 


PEPF B/L-2 


Mass, kilograms 

Diameter, meters 

Drag coefficient, C_ . 

D,D 

Suspension-line unfurling resistance force, newtons 
Canopy unfurling resistance force, newtons 
Mortar ejection velocities, meters/second 


0.9 

0.8 

.3 

.3 

.8 

.8 

8.9 

8.9 

26.7 

26.7 

40.4 

36.6 


Parachute Configurations 

The PEPP B/L-2 flight test parachute was a disk-gap-band configura- 
tion having a nominal diameter of 19*7 meters. The linear -mass-density 
distribution for this parachute was calculated from detailed construc- 
tion diagrams which are given in reference lL. This calculated mass 
distribution, for which it is assumed that the suspension system 

consists of suspension lines only, is shown in figure ll»a. The BLDT ; 

t- 

AV-U flight test parachute was also a disk-gap-band configuration and 

* 

had a nominal diameter of 16.15 meters. The linear-mass-density 
distribution for this parachute was calculated from unpublished 
detailed construction blueprints. Figure lUb shows the calculated 
AV-U distribution. The spikes in the calculated mass distributions 
denote areas in which there are severed overlapped layers of canopy 
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cloth and a layer of fabric tape sewn onto the cloth layers around 
the circumference of that canopy section for reinforcement purposes. 

In the present solution method, these spikes are approximated by 
half-sine curves which attain the appropriate maximum amplitude at the 
mid-points of the spikes. Physical characteristics of the two para- 
chutes which are vised in the present simulation are summarized in 
Table 17. 
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TABLE IV 

PHYSICAL CHARACTERISTICS OF PARACHUTES 


Nominal diameter, meters 

Number of suspension lines, n ^ 

Total mass, kilograms 

Suspension-line length, L , , meters 

si 

Total unfurled length of 
parachute, L^, meters 


PEPP B/L-2 

19.7 
72 

32 . k 

23.8 

3^.3 


cLDT AV-h 
16.15 
18 

36.6 

30.1 

38.8 


Suspension-Line Elastic Characteristics 
Suspension lines for the PEPP B/L-2 parachute were fabricated 
from coreless braided dacron cord having a rated tensile strength 
of 2 U 5 O newtons. Suspension lines for the BLDT AV-4 parachute were 
fabricated from similar cord having a rated tensile strengtn of 3900 
newtons. Force-strain curves for representative samples of each type 
cord were obtained in laboratory tests using an Instron tensile- 
testing machine at a near-zero strain rate. These curves are shown in 
figure 15. A value of suspension-line damping coefficient was 
determined parametrically. As the tension at the vehicle never became 
zero during either of the flight-test deployments being studied, the 
damping coefficient was selected as that value required to maintain 
tension at the vehicle during both of the computed deployment histories. 
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Force, newtons/line 


(a) PEPP B/L-2 Suspension lines 



Figure 15.- Force-strain curves for parachute suspension lines 
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Those particular values were determined, to be O.OU 65 N-sec/line for 
PEPP B/L-2 and O.OT N-sec/line for BLDT AV-U. 
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